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We consider the problem of searching for continuous gravitational wave (CW) sources orbiting 
a companion object. This issue is of particular interest because, the Low Mass X-ray Binaries 
(LMXB's), and among them Sco X-1 - the brightest X-ray source in the sky - might be marginally 
detectable with f» 2 years coherent observation time by the Earth-based laser interferometers ex- 
pected to come on line by 2002, and clearly observable by the second generation of detectors. 
Moreover, several radio pulsars, which could be deemed to be CW sources, are found to orbit a 
companion star or planet, and the LIGO/VIRGO/GEO600 network plans to continuously monitor 
such systems. 

We estimate the computational costs for a search launched over the additional five parameters 
describing generic elliptical orbits (up to e ;S 0.8) using match filtering techniques. These techniques 
provide the optimal signal-to-noise ratio and also a very clear and transparent theoretical framework. 
Since matched filtering will be implemented in the final and the most computationally expensive 
stage of the hierarchical strategies, the theoretical framework provided here can be used to determine 
the computational costs. In order to disentangle the computational burden involved in the orbital 
motion of the CW source, from the other source parameters (position in the sky and spin-down), 
and reduce the complexity of the analysis, we assume that the source is monochromatic (there is no 
intrinsic change in its frequency) and its location in the sky is exactly known. The orbital elements, 
on the other hand, are either assumed to be completely unknown or only partly known. 

We provide ready-to-use analytical expressions for the number of templates required to carry out 
the searches in the astrophysically relevant regions of the parameter space, and how the compu- 
tational cost scales with the ranges of the parameters. We also determine the critical accuracy to 
which a particular parameter must be known, so that no search is needed for it; we provide rigor- 
ous statements, based on the geometrical formulation of data analysis, concerning the size of the 
parameter space so that a particular neutron star is a one-filter target. This result is formulated in 
a completely general form, independent of the particular kind of source, and can be applied to any 
class of signals whose waveform can be accurately predicted. 

We apply our theoretical analysis to Sco X-1 and the 44 neutron stars with binary companions 
which are listed in the most updated version of the radio pulsar catalogue. For up to f» 3 hours of 
coherent integration time, Sco X-1 will need at most a few templates; for a week's integration time 
the number of templates rapidly rises to ~ 5 x lO". This is due to the rather poor measurements 
available today of the projected semi-major axis and the orbital phase of the neutron star. If, 
however, the same search is to be carried out with only a few filters, then more refined measurements 
of the orbital parameters are called for - an improvement of about three orders of magnitude in the 
accuracy is required. Further, we show that the five NS's (radio pulsars) for which the upper-limits 
on the signal strength are highest, require no more than a few templates each and can be targeted 
very cheaply in terms of CPU time. 

Blind searches of the parameter space of orbital elements are, in general, completely un-affordable 
for present or near future dedicated computational resources, when the coherent integration time is 
of the order of the orbital period or longer. For wide binary systems, when the observation covers 
only a fraction of one orbit, the computational burden reduces enormously, and becomes affordable 
for a significant region of the parameter space. 
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I. INTRODUCTION 



The construction of several large-scale interferometric gravitational wave detectors, with optimal sensitivity in the 
frequency window ~ 10 Hz - 1 kHz is close to completion. For a period extending upto a year, starting from the end 
of year 2000, engineering data runs will be carried out in order to test and debug the detector components. Finally, 
between 2002 and 2004 the interferometers will carry out the first set of "science" data-runs with the realistic goal of 
directly observing gravitational waves (GW's). This initial phase will be followed by substantial upgrades on most of 
the instruments aimed at reaching better sensitivity and larger observational bands. The projects include the Laser 
Interferometer Gravitational-wave Observatory (LIGO) in the USA consisting of two facilities, one at Hanford (WA) 
and the other at Livingston (LA), hosting two 4-km and one 2-km interferometers [Q; VIRGO, a French/Italian project 
located at Cascina, near Pisa, Italy, which consists of a 3-km detector (and currently running about one-to-two years 
behind the time-frame presented above) ||^; GEO600, a German/British 0.6-km interferometer located at Ruthe, near 
Hannover, Germany; in Japan, the TAMA project is currently running a medium scale interferometric detector of 
arm length 300 meters, and is planning to extend the baseline to 3 km and carry out other substantial improvements 
on the instrument within the next few years ; finally, the ACIGA consortium will build, if funding is approved, a 
500 meter interferometer (AIGO500) near Perth, Australia [Q. In the meantime, a number of existing resonant bar 
detectors are steadily increasing their sensitivity in narrow bandwidths (~ 1 Hz) covering the kHz spectral window |^ 
and serious efforts are being made to fly, possibly by 2010, a space-borne laser interferometer [Q (LISA: the Laser 
Interferometer Space Antenna). LISA will open the low frequency window 10"*^ Hz - 10^^ Hz, currently accessible 
only via the Doppler tracking of interplanetary space-craft g, and is projected to have sensitivity several orders of 
magnitude better. 

Building and running such powerful machines represents an enormous enterprise, which has started over 30 years 
ago. However, analyzing the large amount of data - several Mbytes per second - and digging out with high confidence 
astrophysical signals from the noise which severely corrupts the data, presents its own challenges. In fact, over the 
past few years, interest has been growing in the area of GW data analysis, which is now regarded as one of the key 
aspects for the successful detection of GWs. 

Several types of GW sources have been envisaged which could be directly observed by Earth-based detectors 
(see [p|-p^ and references therein for recent reviews): (i) Burst sources - such as binary systems of neutron stars 
(NS) and/or black holes (BH) in their in-spiral phase, BH/BH and/or BH/NS mergers, and supernovae explosions 
- whose signals last for a time much shorter, typically between a few milli-seconds and a few minutes, than the 
planned observational time; (ii) stochastic backgrounds of radiation, either of primordial or astrophysical origin, and 
(iii) continuous wave (CW) sources ~ e.g. rapidly rotating neutron stars - where a weak deterministic signal is 
continuously present in the data stream. 

In this paper we deal with data analysis issues regarding CW sources. More specifically, we investigate the com- 
putational load involved in filtering the data stream to search for monochromatic radiation emitted by a neutron star 
orbiting a companion object. 



A. The data analysis challenge 

CW emitters pose one of the most computationally intensive problems in GW data analysis. The weakness of the 
expected signal requires very long observation times, of the order of a year (or possibly more) for accumulating enough 
signal-to-noise ratio (SNR) so that detection is possible. During the observation time, a monochromatic signal in the 
source reference frame, is Doppler modulated by the motion of the detector carried by the spinning Earth orbiting 
the Sun. The emitted energy is spread over ~ 2 x 10^ (T/10''sec)^ (//I kHz) frequency bins of width A/ — l/T, 
where T is the time of observation (the formula holds for T upto six months, after that the number of bins increases 
linearly with T). In order to recover the whole power in one frequency bin one has to "correct" the recorded data 
stream for each possible source position in the sky. The problem is made worse, if the the intrinsic frequency of the 
source changes, say due to spindown. Then the power is spread over 3 x 10^ (r/lO^ yrs)""'^ (T/10^ sec)^ (//I kHz) bins, 
where t = /// is (twice) the spin-down age of the NS. Indeed, one then needs to correct, in addition, for this effect 
searching through one or more spin-down parameters. It is clear that searches for CW's are limited by the available 



computational resources [13 1^. The optimal technique of tracking coherently the signal phase - matched filtering 
- is not affordable for integration times longer than ^ 1 to 10 days (the integration time depends on the range of 
parameters over which the search is being carried out). So far, the only viable strategy for a "blind" search over a 
wide range of parameters is based on a hierarchical structure ||l^,|l^ . In this type of search, coherent and incoherent 
search stages are alternated in order to identify candidate signals - and the relevant parameter range - " cheaply" by 
a suboptimal algorithm. Then, the candidate signals are foUowed-up with a coherent and computationally expensive 
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search over the entire observation time - this step constitutes the bottle neck for processing power - to recover the 
optimal SNR. 

However, due to the large computational burden, the algorithms investigated so far have been restricted to isolated 
NS's. The problem of searching for a CW emitter orbiting a companion object is simply considered computationally 
intractable, as it would require upto five more search parameters which would compound the already enormous 
computational cost. 

To gain insight into the additional Doppler modulation introduced on the signal phase by the source motion around 
a companion, let us consider a NS, of mass tons = 1.4M0, with an orbiting companion object of mass TO2 and period 
P. During the observation time T, the NS velocity changes by Av = (27r/P)'*/'^ (Gtons)^^'^ 9/(1 + 'if'^'^T, where 
q = m2/miMs- The frequency is Doppler shifted by f = ±fAv/c which has to be compared to the frequency 
resolution A/. (Here and after, G denotes the Newton's gravitational constant and c denotes the speed of light or 
actually GW). The maximum observational time T^ax after which one needs to correct for the CW's source orbital 
motion is therefore. 
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Of course, this time scale can vary by orders of magnitude, depending on the NS companion and the orbital period. 
For the relevant astrophysical situations we have: 
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For comparison, we note that the maximum integration time before which the motion of the detector from Earth's 
rotation about its axis, produces a Doppler shift in frequency greater than the frequency resolution b in is Trnax ^ — 
50 (//lkHz)~^/2 min. It is then evident that in the case of a CW source orbiting a companion object, the frequency 
shift is typically more severe and the computational challenge that one faces when searching for the additional 5- 
dimensional parameter vector describing the binary orbit is considerable. 

B. Expected sources and motivation 

The characteristic amplitude of the continuous signal emitted by a triaxial NS, rotating about a principal axis with 
period Pr and frequency fr = 1/Pr reads 

,..„xlO-(^)(^)"(^)", ,1.3) 

where eg is the equ ator ial ellipticity of the star, and we have taken the standard value 10^^ gr cm^ for the NS moment 



of inertia. In Eq. ( |l.3D we assume that the GW energy is released at exactly twice the rotation frequency, so that 
The characteristic amplitude has to be compared with the noise fluctuations. 



1/2 / ^ \ -1/2 



h - 10-25 f y { ^ \ HA) 

'^"-V^-'° Uo-«Hz-ij (^Yo^j (1-4) 

in a resolution bin, where S'„(/) is the noise spectral density of the instrument. A signal is detected if he <L fc/in, 
where the multiplicative factor fc, of order of a few, depends on the number of filters used in the analysis and on 
the (subjective) confidence level of the measurement. The strength of the signal depends crucially on the eg, which 
is a measure of the non-axis-symmetry of a star. Although large uncertainties exist in the maximum value that eg 
can achieve, theoretical studies of the braking strain of NS crusts suggest that ^ 10^^. One can observationally 
constrain eg if the NS rotational period Pr and its first time derivative Pr are known: 




^e<5AxlO-M— i- . (1.5) 
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This yields an upper-limit on the amplitude of the emit ted GW's, by assuming that all energy is lost through 
gravitational wave emission. By combining Eqs. (L3) and ( |l.5| ), we can immediately place the upper- limit 
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The comparison of /ic for some known NS's with /i^ of the laser interferometers is shown in Fig. 0. 




FIG. 1. Upper limits on the GW amplitude from known radio pulsars and Sco X-1. The plot shows the upper-limit on 
he, Eq. ( |l.6| ), for the neutron stars included in the latest version of the published radio pulsar catalogue for which the 
spin-down and the distance are known (isolated pulsars: squares; binary pulsars: bold bullets), and the estimated value of 
he for Sco X-1 (triangle), Eq. ([J). We also show the noise amplitude hn for an observation time T = 10 sec, for the first 
generation detectors (GEO600: solid line; LIGO I: dotted line; VIRGO: dashed line). 



In Table g we also show some of the parameters of the binary radio pulsars which are known so far. It is interesting 
to notice that a bunch of systems are within the reach of the first generation of instruments for coherent observation 
times from 15 months to 4 years, and well within the enhanced configuration sensitivity. It is also worth noticing 
that they represent a fair sample of orbital companions, as they include NS/NS binaries, NS/ white-dwarf binaries, 
NS / main-sea uence stars and NS/planet. However, the actual signal strain could be orders of magnitude below the 
upper-limit ( |l.6|) . 

A far more promising scenario is the accretion of hot material onto the NS surface. Here the induced quadrupole 
moment is directly related to the accretion rate, which can be copious. The gravitational energy reservoir, 
moreover, can be continuously replenished, if persistent accretion occurs. The key idea behind this scenario 
is that gravitational wave radiation can balance the torque due to accretion, and was proposed over 20 years 
ago [|l^,|l9|. However, it has attracted considerable new interest in the past two years and has been fully re- 
vitalized by the launch of the Rossi X-ray Timing Explorer, designed for precision timing of accreting NS's. 
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TABLE I. Neutron stars detected as radio pulsars in binary systems. The table shows the main parameters of the 44 known 
neutron stars detected as radio pulsars orbiting a companion object listed in the Taylor et al. catalogue We list the 

gravitational wave frequency / - assumed to be exactly twice the radio frequency - its first time derivative /, the estimated 
distance D from the solar system and the spin-down age (we do not give all significant digits). The orbital parameters - the 
orbital period P, the projection of the semi-major axis along the line of sight ftp and the eccentricity e - are tabulated with 
their respective errors (in parenthesis). Notice that out of the 44 known systems, seven sources radiate at a frequency < 10 Hz 
and therefore are not observable by the (so far) planned ground-based experiments. 
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0024-7204E 


0021-72E 


565.6 




4.5 




2.256844 (3) 


1.978 (3) 


0.000 (3) 


0024-72041 


0021-721 


573.9 


_ 


4.5 


_ 


0.226 (3) 


0.039 (2) 


0.00 (5) 


0024-7204J 


0021-72J 


952.1 




4.5 




0.120665 (2) 


0.0402 (2) 


0.00 (5) 


0034-0534 




1065.4 


3.8 


1.0 


9.6 


1.58928180 (3) 


1.437768 (5) 


0.0000 (1) 


0045-7319 




2.2 


10.5 


57.0 


6.5 


51.169226 (3) 


174.2540 (8) 


0.807995 (5) 


02184-4232 




860.9 


27.8 


5.8 


8.7 


2.02885 (1) 


1.9844 (1) 


0.00000 (2) 


0437-4715 




347.4 


3.4 


0.1 


9.2 


5.741042329 (12) 


3.3666787 (14) 


0.0000187 (10) 


0613-0200 




653.2 


2.3 


2.2 


9.6 


1.19851251 (2) 


1.09145 (2) 


0.000000 (22) 


0700+6418 


0655+64 


10.2 


0.04 


0.5 


9.7 


1.028669703 (1) 


4.125612 (5) 


0.0000075 (11) 


0751+1807 




574.9 


1.3 


2.0 


9.8 


0.2631442679 (4) 


0.396615 (2) 


0.0000 (1) 


0823+0159 


0820+02 


2.3 


0.3 


1.4 


8.1 


1232.47 (12) 


162.1466 (7) 


0.011868 (7) 


1012+5307 




380.5 


1.1 


0.5 


9.8 


0.604672713 (5) 


0.581816 (6) 


0.00000 (2) 


1022+10 




121.6 




0.6 




7.80513014 (8) 


16.765412 (6) 


0.0000977 (6) 


1045-4509 




267.6 


0.7 


3.2 


9.8 


4.0835291 (1) 


3.015107 (9) 


0.000019 (6) 


1300+1240 


1257+12 


321.6 


5.9 


0.6 


8.9 


66.536 (1) 


0.0013106 (6) 


0.0182 (9) 


1302-6350 


1259-63 


41.9 


1994.5 


4.6 


5.5 


1236.72359 (5) 


1296.580 (2) 


0.869931 (1) 


1312+1810 


1310+18 


60.3 




18.9 




255.8 (6) 


84.2 (7) 


0.002 (3) 


1455-3330 




250.4 


_ 


0.7 


_ 


76.17458 (2) 


32.36222 (6) 


0.000167 (2) 


1518+4904 




48.9 


_ 


0.7 


_ 


8.6340047 (8) 


20.043982 (8) 


0.2494848 (7) 


1537+1155 


1534+12 


52.8 


3.4 


0.7 


8.4 


0.4207372998 (3) 


3.729468 (9) 


0.2736779 (6) 


1623-2631 


1620-26 


180.6 


12.9 


1.8 


8.3 


191.442819 (2) 


64.809476 (13) 


0.0253147 (5) 


1640+2224 




632.2 


0.6 


1.2 


10.2 


175.460668 (2) 


55.329722 (2) 


0.00079742 (5) 


1641+3627B 


1639+36B 


566.9 




7.7 




1.259113 (3) 


1.389 (2) 


0.005 (5) 


1643-1224 




432.7 


3.1 


4.9 


9.3 


147.01743 (6) 


25.07260 (4) 


0.000506 (2) 


1713+0747 




437.6 


0.8 


1.1 


9.9 


67.82512988 (2) 


32.342413 (1) 


0.00007492 (2) 


1721-1936 


1718-19 


2.0 


3.2 


11.6 


7.0 


0.2582735 (2) 


0.3526 (8) 


0.000 (5) 


1748-2446A 


1744-24A 


173.0 


-0.3 


7.1 




0.0756461170 (4) 


0.11963 (2) 


0.0000 (12) 


1803-2712 


1800-27 


6.0 


0.3 


3.6 


8.5 


406.781 (2) 


58.9397 (18) 


0.000507 (6) 


1804-0735 


1802-07 


86.6 


1.8 


3.2 


8.9 


2.6167634 (5) 


3.92047 (4) 


0.211999 (15) 


1804-2718 




214.1 




1.2 




11.15 (1) 


7.20 (4) 


0.00 (1) 


1823-1115 


1820-11 


7.1 


35.2 


6.3 


6.5 


357.7622 (3) 


200.672 (3) 


0.79462 (1) 


1834-0010 


1831-00 


3.8 


0.1 


2.6 


8.8 


1.811103 (2) 


0.7231 (8) 


0.000 (4) 


1857+0943 


1855+09 


373.0 


1.2 


0.9 


9.7 


12.3271711905 (6) 


9.2307802 (4) 


0.00002168 (5) 


1910+0004 




552.7 




4.1 




0.140996 (1) 


0.038 (2) 


0.00 (1) 


1915+1606 


1913+16 


33.9 


4.9 


7.1 


8.0 


0.322997462736 (7) 


2.3417592 (19) 


0.6171308 (4) 


1955+2908 


1953+29 


326.1 


1.6 


5.4 


9.5 


117.349097 (3) 


31.412686 (5) 


0.0003304 (3) 


1959+2048 


1957+20 


1244.2 


13.0 


1.5 


9.2 


0.3819666389 (13) 


0.892268 (8) 


0.00000 (4) 


2019+2425 




508.3 


0.9 


0.9 


9.9 


76.5116347 (1) 


38.767625 (1) 


0.00011111 (6) 


2033+17 




336.2 




1.4 




56.2 (1) 


20.07 (8) 


0.00 (5) 


2130+1210C 


2127+llC 


65.5 


10.7 


10.0 


8.0 


0.335282052 (6) 


2.520 (3) 


0.68141 (2) 


2145-0750 




124.6 




0.5 




6.83890256 (8) 


10.16411 (1) 


0.000021 (2) 


2229+2643 




671.6 


0.4 


1.4 


10.4 


93.015895 (2) 


18.912519 (5) 


0.0002550 (4) 


2305+4707 


2303+46 


1.9 


1.0 


4.4 


7.5 


12.33954454 (17) 


32.6878 (3) 


0.658369 (9) 


2317+1439 




580.5 


0.4 


1.9 


10.4 


2.459331464 (2) 


2.3139483 (9) 


0.0000005 (7) 



The observational evidence that Low Mass X-ray Binaries (LMXB's) - binary systems where a compact object ac- 
cretes material from a low mass companion - in our Galaxy are clustered around a rotation frequency ~ 300 Hz, led 
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Bildsten [ po| to propose a mechanism to explain this behaviour. The fundamental idea is that continuous emission of 
GW's radiates away the angular momentum that is transfered to the NS by the infalling material. The fact that the 
rate of angular momentum loss through GW's scales as provides a very natural justification of the clustering of 
rotation frequency of several sources. The physical process responsible for producing a net quadrupole moment is the 
change of composition in the NS crust, which in turn is produced by the temperature gradient caused by the in- falling 
hot material. Recently, Ushomirsky et al. have posed this initial idea on more solid theoretical grounds. If such 
mechanism does operate, LMXB's are extremely interesting candidate sources for Earth-based detectors. Several 
systems would be detectable by LIGO operating in the "enhanced" configuration (LIGO II), if the detector sensi- 
tivity is tuned, through narrow-banding, around the emission frequency. In particular, Sco X-1, the most luminous 
X-ray source in the sky, possibly is marginally detectable by "initial" LIGO and GEO600 (the latter in narrow-band 
configuration), where an integration time of approximately 2 years would be required. The characteristic amplitude 
for this class of sources is: 

/i,~4x 10-27 (10-8 erg cm-2 sec- ij ' ^^'^^ 

where Rns is the typical radius of a neutron star and Fx is the X-ray flux. 

Data analysis strategies to search for NS's in binary systems might therefore be vital for ground-based detectors. 
They would supplement the already existing strategies for the several sources mentioned above. 



C. Organization of the paper 



The paper is organized as follows: In Section || we set up the parametrized model of the signal based on the two 
body problem. We then identify the search parameters of the signal, which are essentially the binary orbital elements, 
for the general case of elliptic or bits, a s well as for circular orbits. In Section III we discuss the geometrical approach 
to the data analysis; in Section III A we review the method of determining the number of templates which involves 



the computation of the metric on the si gnal manifold, and we also state useful properties of the metric which are then 



used in the later sections. In Section HI B we present a rigorous criterion, naturally suggested by the geometrical 
picture, for estimating the number of filters required in targeted searches and to determine whether any particular 
source is a one-filter target. 

Sections IV and contain the key results of the paper. In Section IV we obtain analytic approximate formulae 
for the number density of filters per unit coordinate volume in the parameter space in the two important regimes of 
observation time: (i) observations lasting for several orbital periods, and (ii) observations lasting a fraction of an orbit. 
These expressions provide for us the scaling: that is, how the computational cost scales with different parameters. 
They can also be readily applied to estimate the computational cost of detecting GW from some observed systems. 
Further, we also compare these results with the fully numerical computation of the metric determinant (which shows 
excellent agreement over a wide range of the parameters) , and determine the range of applicability of these approximate 
expressions as a function of the relevant parameters (in particular the time of observation as compared to the orbital 
period). In the Section ^ we apply the above analysis to known NS's in binary systems, specifically radio-pulsars and 
LMXB's. Given the current (or near future) computational resources, we then determine the systems that could be 
searched by matched filtering methods. Section VI contains our final remarks and future directions of work. 



II. THE MODEL 



A. The description of the problem 



In this paper we employ the matched filtering technique for coherent integration of the signal. It provides for us the 
optimal signal-to-noise ratio that can be achieved with linear data analysis techniques, and also an elegant theoretical 
framework can be given in which the computational costs can be derived in a very transparent way. We have already 
remarked that the processing power of present /near future computers falls woefully short for carrying out full matched- 
filter based searches of CW sources. However, it is important to note that matched filters are implemented in the 
final stage of every hierarchical strategy that has been devised so far [ p^p6| . Our analysis can be therefore applied 
to the typically reduced parameter range which is relevant in this final step. Since the final stage accounts for most 
of the processing power required in the entire search, our results essentially reflect the computational costs involved 
in the search. 
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The amplitude of the GW's is expected to be small for CW sources. Therefore long observation times are imminent 
for extracting the signal out of the detector noise with a reasonable signal-to- noise ratio (SNR). The observation times 
could last for a few months, up to a year or more. During this period of observation it becomes essential to take into 
account the changes in frequency which result from the variation in relative positions of the source and the detector 
- Doppler effects ~ and also those that arc intrinsic to the source, such as the spin-down. The general idea is (i) to 
correct for these changes in the phase, (ii) use the FFT algorithm to compute efficiently the power spectrum and 
then, (iii) look for a statistically significant peak in the power spectrum itself. 

The aim of this paper is to estimate the extra computational costs associated with the search for sources in binary 
systems. In order to disentangle the phase Doppler modulation produced by the source motion around the orbital 
companion, and limit the complexity of our analysis, we will make the following assumptions: 

• The GW source radiates perfectly monochromatic GW's at the intrinsic frequency /o, that is the frequency in 
its own rest frame is constant: we ignore the effect of spin-down parameters. 

• The location of the source in the sky is perfectly known: the Doppler effect produced by the motion of the Earth 
spinning on its own axis, orbiting the Sun etc. can be computed and subtracted from the data stream with very 
little CPU load. 

These hypotheses leave us only to grapple with the problem of the motion of the source. Notice that they do not 
affect the generality of our analysis. This approach returns the extra cost involved in searching for NS's in binary 
systems, for each template of the isolated NS search (with the appropriate number of spin-down parameters and 
range of parameter values). Our analysis thus provides a solid upper-limit on the total computational costs (it ignores 
possible correlations between the parameters). In fact, if we call A/is the number of templates to search for an isolated 
NS (the search parameters are the location in the sky and spin-down coefficients), and Af the number of templates 
required to search for a perfectly monochromatic source in binary orbit whose location in the sky is assumed to be 
known, the search parameters being only the orbital elements of the binary system, then the total number of filters 
A/t required for a generic source, is bound by the following inequality: 

AAt < Ms X A/". (2.1) 

We take the orbit of the source to be Keplerian, elliptical in shape, but we assume that we do not exactly know, or 
do not know at all, the orbital parameters. It is possible that a particular parameter or a set of parameters is known 
so accurately that no search is required over these parameters. Such a situation arises for some of the known radio 
binary pulsars listed in the published catalogue We then have a targeted search, where one must now search 

only over rest of the parameters. In the extreme case, all parameters may be known so accurately, that only a single 
Doppler correction needs to be applied (one filter in the language of data analysis). Usually for known systems, some 
orbital elements are known within given error bars. For example, one might know the period of the orbit within 
certain limits, then these limits provide for us the ranges of the parameters, over which the search must be carried 
out. On the other hand we might not know anything at all about the binary system, in which case one must launch 
a more extensive search with a large number of templates covering the entire region of the parameter space - this is 
the so-called "blind" search. 



B. The two body problem, orbital elements and search parameters 

The essential problem of celestial mechanics is the two-body problem, where one solves for the motion of two point 
particles of masses mi and m2, and position vectors ri and r2, attracted by their mutual gravitational force. The 
solution is simple in their centre-of-mass frame in terms of the relative position vector r = ri — r2. The vector r lies 
in a plane and traces out an ellipse; the individual masses trace out similar ellipses scaled by factors depending on 
the masses: 

ri = r, ^2 — r, (2.2) 

where M = nii + m2- Here, we (arbitrarily) take ri = to describe the position of the CW source we want to 
detect. We assume that the centre-of-mass of the binary system can be considered at rest, or in uniform motion, with 
respect to the Solar System barycentre during the period of observation; therefore knowing ri as a function of time 
is sufficient for us to compute the GW Doppler phase shift (f>D{t)- 

It is convenient to set up a Cartesian coordinate system (<^,77, C) attached to the binary source, so that: (i) ri lies 
in the (^, rj) plane with the origin at the centre of the ellipse; (ii) the semi-major axis of the ellipse coincides with 
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the ^-axis, and (iii) the (^-axis is perpendicular to the orbital plane and points in the direction of the orbital angular 
momentum. We specify the direction to the detector, in the 77, C) frame, by the unit vector: 



n = (sin e cos tjj) ^ + (sin e sin -0) 17 + (cos e) C ■ 



(2.3) 



where e and ■0 are the usual polar angles. 

The Doppler phase shift at the detector output therefore reads: 

0D(t) = -^ [fi-ri(t)] , 



(2.4) 



where c ~ 2.9979 x 10^° cm/sec is the speed of light. 

The orbit in the (^, 77) plane is given as follows: Let a be the semi-major axes of the elliptical orbit and e the 
eccentricity, then, the orbit is described by the equations: 



i{t) = a cos E{t), 

r](t) = ayj\ ~ e^ smE{t), 



(2.5) 



where E is the so-called eccentric anomaly, and is a function of the time t. It is related to the mean angular velocity 
Lu and the mean anomaly M{t) by the Kepler equation. 



E{t) - e sin E{t) = ujt + a = M{t), 



(2.6) 



where a is an initial phase, Q < a < 2-k. When ut + a — Q we have E — Q and the mass is closest to the focus 
^ = ae,?7 = 0. These equations describe the orbit in the (^,7y)-plane as function of time, that is determined by the 
four orbital elements a,ijj,a,e. However, the orbit in space requires two additional parameters, namely, the angles e 
and V' mentioned before. Thus, in all we have six orbital elements which specify the orbit in space. 
The Doppler phase correction, Eq. (p]j), is obtained from Eq. (2.5), and reads as follows: 



27r/oasine 



iipcosE{t) + smijj\/l — e^ s\nE{t) 



(2.7) 



The total phase <I> at the barycentre (we have assumed that the corrections have been made for the Earth's motion) 
is therefore, 



$(t) = 27r/oi + 0i5(t). 



(2.8) 



Since we have assumed that we have corrected for the Earth's motion, the t can be regarded as the barycentric time. 

For computing the metric in the next section we will require the time averages of the derivatives of the phase with 
respect to the search parameters. Since the Kepler equation connects E to the time t, or M implicitly, cosE and sin_B 
are usually expressed as power series in the eccentricity parameter e with harmonics in M . Also the time averages 
are more conveniently computed using M. We thus write formally, 



cos E = ^Ck{e) cos(fcM) , 



k=0 



\/l-e2 sinE Sk (e) sin(fcM) , 



(2.9) 



fe=i 



where Cfc(e) and Sk{e) are power series in e. In this work we consider expansions up to the 7th power in e, which 
automatically means that we consider 7 harmonics in M. The Ck{e) and Sk{e) are given in Appendix ^ upto this 
order. 

The parameters over which one must launch a search are not exactly the orbital elements, and need not be of the 
same number. It is the Doppler phase correction that is observed and so the information about the system that we 
can glean depends on the combination of the orbital elements that enter into it. The a and e combine into a single 
parameter a sine = Op, the projected semi-major axis along the line of sight, which is actually the quantity inferred 
from astronomical observations. The other search parameters are the remaining orbital elements uj, a, e and ip. So in 
the general case, when we do not know any of the parameters exactly, we have a 5-dimensional parameter space to 
search. 
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We note that fo is not a search parameter because of the special search technique that is employed |14|: it involves 
the "stretching" of the time coordinate in such a way, so as to make the signal appear monochromatic in this time 
coordinate. One then simply takes the FFT to compute the power spectrum, which now is concentrated in a single 
frequency bin. 

The gravitational waveform that we are dealing with, reads (in the barycentric frame): 

hit, A) = 5R[ylexp(-i$(i; A) + i*)], (2.10) 

where $(t; A) is given by Eq. ( ^.8[ ), A — {ap,uj, a, e, "0) is the 5-dimensional vector that refers to the parameters that 
are required in the discrete mesh search, and A = (/o. A). The polarization amplitude A, and the polarization phase 
^' are slowly varying time-dependent functions over the time-scale of the day and depend on the relative orientations 
of the source and the detector. In agreement with all the investigations carried out so far, we assume them constant 
in our analysis. It is expected that these factors can be easily included in the full analysis and will not significantly 
affect the computational burden Il4]. 



III. THE DATA ANALYSIS 

In the geometrical picture p3| , p^ , the signal is a vector in the vector space of data trains and the iV-parameter 
family of signals traces out an A'^-dimcnsional manifold which is termed as the signal manifold. The parameters 
themselves are coordinates on this manifold. One can introduce a metric jjk on the signal manifold which is related 
to the fractional loss in the SNR when there is a mismatch of parameters between the signal and the filter. The 
spacing of the grid of filters is decided by the fractional loss due to the imperfect match that can be tolerated. Given 
the parameter space that one needs to scan, it is then easy to estimate the total number of filters required to carry out 



the search for the signal. In Sec. Ill A we first briefly review the method introduced by Brady et al. p4f| which in turn 



was based on Owen's [E6[ method for searching for GW signals from in-spiralling compact binaries. We also pres ent 



we 



some useful general properties of the metric on the signal manifold which are used in the next section. In Sec. [II B 
then consider the case where accurate information about one or more parameters of a NS are available (for instance 
through radio observations) so that the parameter space that one needs to search over, is drastically reduced. We 
present a rigorous approach, based on differential geometric methods, to determine the exact dimensionality of the 
data analysis problem and formulate criteria to decide whether a given source is a one-filter target. 



A. Number of filters: General 



In the method presented in the key idea is to first correct for the Doppler effect in the phase of the signal for 
each of the grid points of the parameter space and then compute the power spectrum. The latter is obtained efficiently 
via the FFT algorithm. Even if we know the frequency of the pulsar, it is desirable to search over a band ~ 1% of 
the pulsar frequency pst . We therefore have a large number of frequency bins to search over and the FFT algorithm 
is thus computationally advantageous. If the Doppler correction is right, that is, if the signal and filter parameters 
match perfectly, then the signal is all concentrated at / = /o in the power spectrum. The grid spacing is decided 
by the amount the maximum of the power spectrum falls, when the parameters of the signal and filter mismatch. 
The mismatch fj, is defined as the fractional reduction in the maximum of the power spectrum when the parameters 
mismatch. Fixing the mismatch ^, fixes the grid spacing of the filters in the parameter space which we will denote by 
V. The number density of filters (the number of filters per unit proper volume - proper volume defined through the 
metric) in V depends on /i, and is denoted by p^ifJ'), where N is the dimension of V. For a hyper-rectangular mesh 
(which does not represent necessarily the most efficient tiling of the parameter space) it reads: 




In fact, the proper distance dl between two filters is: 




The proper volume of V, can be easily computed from [ p^ : 
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Vv = J^dX ^det||7,fe||; 



the number of filters Af is then just the proper volume (3.3), times the filter density (3.1) 



(3.3) 
(3.4) 



For the elliptical orbit we have N = 5, while for the circular case the number reduces to = 3. See Table ^ for the 
relevant values of pjv(/^)- 

TABLE 11. The number density of filters pjv(At), Eq. (|t|) for the circular and elliptical orbit case (the dimensions of the 
parameter space are A*' = 3 and N = 5, respectively) as a function of the mismatch fi. 



dimension of the 


mismatch 


parameter space 


1% 


3% 


5% 


10% 


30% 


3 


649.5 


125.0 


58.1 


20.5 


3.9 


5 


174963.0 


11206.5 


3125.0 


552.4 


35.4 



If the signal parameters are A — (/o, A) and the filter parameters are A + AA, the power spectrum for an observation 
time T is given by, 



P(/;/o,A,AA) 



T 



dt exp[i$(t; A, AA)] 



where 



$(i; A, AA) - 27r(/ - fo)t + cf^oit; /o, A + AA) - ^oit; fo, A) 



(3.5) 



(3.6) 



The mismatch is both in A as well as in / (this can occur because of sampling at the wrong frequency) and is denoted 
by m(A, AA): 



m(A, AA) ^ 1 - ffc^MA) AA" AA/^ + o{AA') . 

^(JO, Jo, A, Uj 



From Eqs. (3^) and (3/7) the metric gap can be computed by Taylor expansion. It is given by, 

= ($c$/3) - (*a)(*/3), 



(3.7) 



(3.8) 



where the suffix, say a, denotes derivative with respect to AA" and the angular brackets denote time averages defined 
as follows: For a function X{t) defined on the data train [0,T], the time average of X is. 



1 
T 



dtX{t). 



(3.9) 



We remark that gap is not the metric which is used to calculate the proper volume, because it still includes /. We 
need to maximize over /, which is tantamount to projecting gaf3 orthogonal to the A/ direction. Thus the metric on 
the sub-manifold of the search parameters A is. 



7jfe = gjk 



9Qj90k 

500 



(3.10) 



Here the Greek and Latin indices range over the parameter A (a, /3 = 0, 1, ■■■,N) and A (j, k — 1,2, N), respectively; 
the index identifies the parameter corresponding to the frequency /. If /max is the highest GW frequency that we 
are searching for, then we must put fo = /max in the above expression f or 7 ,fc. T he p roper volum e of the signal 
and the total number of templates can be easily derived by inserting Eqs. (3T) and ( 3.10 ), into Eqs. ( |3.3| ) and ([3.4|), 
respectively. Notice that the parameter fj, is the projected mismatch, where the maximization over / has already been 
accounted for. 

The metric defined above through the phase <E> has certain elegant properties which we will use in the next section 
to simplify computations: 
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1. Scaling: If we scale $ by a constant factor that is if we define, 

$ = X^, (3.11) 
then each component of the metric ga0 and ^jk is scaled by the factor x^' 

gap = gaf} , ijk = lok ; (3.12) 

the determinants are scaled accordingly: 

det||5a/3|| ^x'^^^' det||ff„;3||, det ||7,fe|| = x'"^ det |fe|| ; (3.13) 
the proper volume V-p is scaled by the square root of the determinant, and so is M: 

Vv = {x'')Vv, AA-(x^)AA. (3.14) 

2. Translation: If a function only of the parameters, say /(A), not containing time is added to the metric 
remains invariant under this transformation. 

These properties are easily verified. We use the first property to make all the coordinates dimensionless in the metric, 

compute Vp, the proper volume in the scaled metric, and then finally multiply by the appropriate power of x to 
obtain the actual volume V-p and also the number of filters. We have already used the second property to shift the 



origin from the focus to the centre of the elliptical orbit in Eq. (2_^). The expressions for the phases differ by an 
amount proportional to ae, which leaves the metric unaltered. 

B. Targeted searches 

So far we have assumed that the size of the parameter space V along any direction is much larger than the distance 
between two filters; as a consequence A*" ^ 1. This is the usual situation which applies when the source parameters 
are not known in advance. However, about a thousand NS's are known today - mainly in the radio band - and in 
looking for GW's emitted by such objects the size of V is drastically reduced, as some or all of the parameters are 
known a priori with a fair degree of accuracy. 

We define: 

ALn<A^<A^„_ J-l,..,iV (3.15) 
the parameter range over which one needs to carry out a search, so that the "error bar" is: 

AA-'" = A^„a, - ALn j = l,2,..,iV. (3.16) 

We distinguish then two cases: 



Blind searches: here A^^x ^ -^minj ^"^^ AA-' ~ A^j^x! -^max ^"^^ ^Iniw ^'^^ decided by the observer, based on the 
available theoretical understanding of the astrophysical scenario; 



• Targeted searches: A^ax ~ '^min' ^"^^ AA^A-' < 1, or even ^ 1; A^^ax 3'iid X^^^^^ are obtained from electromagnetic 
observations of the source. 

In targeted searches, the analysis presented in the previous section for the estimation of the number of parameters 
cannot be applied directly. It must be preceded by the determination of the actual search parameters - which ones 
and how many ~ in other words the number of dimensions of the signal manifold. In fact, some of the parameters, 
might be known so accurately - where the required accuracy is determined by the loss of SNR that cannot exceed 
the maximum mismatch - that one does not need to search through these parameters at all. Of particular interest 
is to determine whether any known NS is a one- filter target, that is A/" = 1. If this condition is satisfied, we need just 
one filter, constructed with the parameter values provided by the available observations. 



It is easy to provide an intuitive example that shows the danger of applying the analysis of Sec. [II A, without first 
checking the effective number of dimensions of the signal manifold. Assume, for sak e of s implic ity, t hat a generic GW 
signal h depends only on two parameters, A^ and A^, with errors given by Eqs. ( 3.15| ) and ( 3.16 ); the metric 7jfc, 
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assoc iate d with this particular waveform, is represented by a 2 x 2 matrix. The number of filters given by Eqs. ( |3.3| ) 
and (3^) is A/'(AA^, AA^), where we have explicitly included the size of the parameter space in the argument of N . 
It is clear that by making, say, AA^ small enough, we can always obtain M < 1. In other words our uncertainty on 
A^ can be arbitrary large, but it seems that by refining our knowledge on A^ we can search for that signal with only 
one filter. This conclusion is clearly wrong; the correct physical interpretation is that, A^ is not a search parameter, 
and we can assign its best value, provided by the available observations, to the entire bank of filters. However, we 
still need to search over A^, which will give us the requisite bank of filters. The key point is therefore to decide what 
is the effective number of dimensions - how many parameters one needs to search for - and then apply the analysis 

of the previous section to search through the appropriate paramet er spa ce. 

Given any signal model h{t] A) and the errors in the parameters ( 3.15| ), ( ^.16| ), we provide here a criterion to decide 
the effecti ve nu mber of dimensions. We consider the generic N-dimensional parameter space described by the metric 
7jfc, Eq. ( 3.1Cl| ), where j^k — 1,2, ...,7V. The proper distance dl between two filters is given by Eq. (^.2|). The key 
quantity is the " thickness'^ of the parameter space along any particular direction A^ associated with AA^ . From 
standard geometrical analysis we can derive it in a straightforward way: 



AA^ 



(3.17) 



Here is the inverse of ^jk, which is related to 7-''^ by 7-''^ — "f^^ /x^ 



is the length of the unit vector orthogonal 
to the hypersurface A-' ~ const. The condition such that a particular A-' is not a search parameter is 



which, via Eqs. (U) and (|t^, becomes 



max [T^] < dl , 



(3.18) 



AA-' < 2 J max 



(3.19) 



If the previous inequality yields for all the N parameters, the source needs at most a few filters. More rigorously, for a 
source to be exactly a one-filter target, the diameter of the parameter space, namely, the maximum distance between 
any pair of points in V, should be less than dl. The parameter space is now a A^-dimensional parallelepiped with 2^ 
vertices. This yields the following condition on the parameter errors: 



(3.20) 



where the maximum is taken over diametrically opposite 2^^^ pairs of vertices of the parallelepiped. If Eq. (3.2C) 
is satisfied, M — 1 and the template is constructed by setting the parameters equal to their best fit provided by the 
observations at hand. In Eq. (3.20) we have tacitly assumed that the size of the parallelepiped is smaller than the 
scale on which the manifold curves. 

If M (< N) parameters satisfy the condition (3.18|), or, equivalently, (3.19), we call ^jk the {N — M) x [N — M ) 



metric that describes the problem, which is constructed using the {N — M) parameters that do not satisfy Eq. ( 3.19| ) 
in the very same way described in the previous secti on. The s igna l manifold is of lower dimension, and the total 
number of filters required can be computed using Eqs. (3.3) and (3.4) for the relevant number of dimensions (N — M): 



dVr 



Vn- 



det I |7jfe I 



(3.21) 



Here the index {N — M) in the density, the number of filters, parameter space, and volume, stresses the fact that we 
are considering a. N — M dimensional manifold; we also use ~ instead of "equal", because of possible edge effects, 
introduced by the dimensions which are close to, but not exactly, known. In the following we discuss more in detail 
this issue, and the subtleties involved in the determination of the number of parameters. 

The picture becomes complicated in the case where the thickness of the parameter space in one or more dimensions 
is smaller than, but of the same order of, the proper distance between two filters: ^ dl. In this case it is useful to 
introduce the following notation: 



= 6dl 



< (5 < 1 . 



(3.22) 
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When 5 ^ 1, we must be able to reach the signal on the boundary of the parameter space, which means that the 
filters should be placed closer than in the case of J <C 1. In the general case, the correct determination of the number 
of filters can be very complicated, due to th e edg e effects. For sake of simplicity, let us consider, first, the case where 
only one parameter, say j — 1, satisfies Eq. ( |3.22| ). In this situation, we have. 



Solving the above equation we have. 



(3.23) 



dl 




(3.24) 



and since the unit of volume in parameter space is dlpatch = dl , the total number of filters is given by 



A/" = pN-iilJ-) 



dVN- 



Vn- 




det ||7jfe| 



(3.25) 



wher e 7, k — 2, .., TV. It is easy to check that we recover the expected formulae for the parameter space volume from 
Eq. ( |3.25 ) in the limiting cases when 5 = 0, 1. 

In more general cases, when two or more parameters are characterized hy 5 ^ 1, the situation becomes very complex, 
and an " exact" determination of the number of filters particularly hard. In the case when the signal manifold is 3- 
dimensional and two parameters, say Ai and A2, are known with 5 1 it is still possible to obtain a rather simple 
expression for the thickness. In this case: 



r 




1/2 



733 



AA^AA* 



(3.26) 



A comparison analogous to the previous one then applies. 

It is however possible to obtain an estimate of the total number of filters which is correct within a factor of a few 
in the general case (we will use this criterion in Sec. ^ when we estimate the computational costs for k nown radio 
pulsars). We introduce the following notation: M (< N) is the n umbe r of param eters that satisfy Eq. ( |3.22 ), and 
M'{< M) those for which 5 <c; 1; the M' parameters satisfy Eq. ( ^.Ig ) or ( |3.19 ), and in practice we can consider 
them as " exactly" known; the value of 5 that we adopt to distinguish the M' parameters from the other M — M' ones 
is 5 = 0.1. Indeed, M — M' is the number of parameters for 0.1 < (5 < 1, those for which edge effects may become 
important. The reason of the choice 5 — 0.1 is the following: in the case M' = M = 1, it is easy to verify, through 
Eq. (3.25), that the estimate of the total number of filters is very accurate. We assume that similar accuracy will 
hold also for few number of dimensions. The values N , M and M' satisfy therefore the inequality Af < M < N. 
To summarize, we consider the M' parameters to be exactly known, M — M' "questionable" (border effects need to 
be taken into account) and N — M to search for. One can imagine the questionable parameters forming a M — M' 
dimensional parallelepiped at each point of the N — M dimensional sub-manifold. The projection of this parallelepiped 
'orthogonal' to the N — M dimensional sub-manifold of search parameters decides the number of filters needed for 
the questionable parameters. The total number of filters is therefore: 



(3.27) 



where M(n-m) is given by Eq. (3.21). The multiplicative factor 2^*^ ^ accounts for the number of filters (which 
we have overestimated) to take care of each of the questionable parameters. Here, we have been cautious to assume 
two filters per parameter, say one corresponding to A^^^jj^ and one to A^^^^. It is simple to check that we overestimate 
the number of filters by applying Eq. ( |3.25| ), which is valid if M' = and M 1. In the case M - M' > 1 Eq. ( |3.27| ) 
might greatly overestim ate th e actual number of filters. In practical cases M — AI' is expected to be at most a few, 
so that the expression (3.27) should provide a meaningful upper-limit (clearly A/at-m < M < 2^^^'^'^ ^J\f[M-M))\ 
however, one can easily resolve this issue by setting up a Monte-Carlo experiment to determine the additional number 
of filters and their location in V for the questionable dimensions. 
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IV. THE COMPUTATIONAL COSTS 



A. The signal phase and dimensionless parameters 



It is convenient for the purposes of computation to express th e ph ase Eg. (2.S), in terms of dimensionless 

parameters and use the above mentioned scaling properties, Eqs. ( 3.11 ) to ( 3.14 ), to obtain the actual proper volume 
Vp, Eq. (O). We write 



$ = (27r/oT) $, 



where 



KU - 



X coaE + Yy/l - e"^ sinE', 



and the dimensionless parameters are given by: 

t 

T ' 

f-fo 
n=——, 

JO 

X = —a sin e cos ip/cT , 
Y = — a sine sin V'/cT, 
VL = ujT . 

Here u is a dimensionless time satisfying < u < 1 . So no w the new set of parameters is A = (k , X , 1" , e , f2 
\= {X ,Y ,e,n,a). The scaling factor x in Eqs. ( |3.12| ), ( |3.13D and (|3.14D is clearly: 



X = 27r/,„axT , 



(4.1) 
(4.2) 

(4.3) 

(4.4) 

(4.5) 
(4.6) 
(4.7) 

a) and 
(4.8) 



wher e w e have set /o — /max, the maximum source frequency that one searches for. Using an equation analogous to 
Eq. (3.8) we obtain the scaled metric gap- 



The exact expression of the determinant of the metric from which we compute the volume is quite complicated. 
However, in this section our goal is to obtain approximate analytical expressions - which turn out to be very accurate 
over most of the relevant parameter range - which give us the important information about how the number of filters 
scales as a function of the key observables. For this purpose it is necessary to make some assumptions about the 
ranges of these parameters from astrophysical scenarios. The semi-major axis ons of a NS orbiting a companion of 
mass 7712 ~ with mass ratio q = 77i2/77iNS - is: 



flNS - 2.65 X 10^ 



"7ns 



1.4 M, 



1/3 



lO-^rad/sec, 



^-2/3 


q 




Ul + g)2/3j 



(4.9) 



where the choice of oj corresponds to ~ 17 hours orbital period (notice that for a NS in a 2 hours orbit, the semi-major 
axis would be a factor ~ 4.2 smaller). We take the time of observation to be typically 4 months long, so that T = 10^ 
sec and we set, as reference, /max = 1 kHz. Thus the factor fmaxT ~ 10^". The dimensionless parameters X and Y 
are X ~ F ~ ons/cT, and give essentially the size of the orbit in units of the distance light (or actually GW) travels 
during the observational time T . They are related to by: 



cT 



8.8 X 10" 



1.4 Mq 



1/3 



17 

103 



-2/3 



T 
107 sec 



-1/3 r 



(1 + g)2/3 



(4.10) 



n^iQ^i / , \ ( -L\ . (4.11) 

VlO-^rad/secy sec ) ^ ' 

Notice that the value of corresponds to an observation that extends over several orbits of the binary system. In fact 
is the number of radians, or r2/27r is the number of orbits, that the star completes in the time T. In this case, one 
can Taylor expand the matrix elements, and therefore the determinant, as a function of <g; 1. Note, moreover, 
that in this case X ^ Y ^ 17^^ <C 1, that is a and w become of the same order - in geometrical units, c = G = 1 
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- only for a NS orbiting a companion object at a distance of the order of the gravitational radius. We will use these 
results in the approximation scheme, retaining only the leading order terms in X and Y . 

The other physically relevant case is in the opposite limit f2 <C 1. We therefore Taylor expand the matrix determinant 
as a function of SI. This case applies to observations that cover only a fraction of one orbital period, and the closed 
form expression that we obtain turns out to provide a rather good approximation of the full expression upto ~ 1. 

In the next subsection we examine the case when the orbit is circular, that is, e = 0. We treat this simpler case 
first in order to obtain useful insights, before we address the general case of the elliptical orbit. 



B. Circular orbits 



The circular case is important both from the pedagogical and physical point of view: (a) it provides us several 
insights into the problem via a comparatively easier computation; (b) in targeted searches, several known NS's in 
binary systems, including Sco X-1, are essentially in a circular orbit; (c) for blind searches, present and near future 
processing power is likely to allow us to search over a reasonable parameter space mainly for emitters orbiting a 
companion with e = 0. 

For circular orbits, the expression of the phase (4.2) simplifies considerably: here e = 0, and ■0 and a combine 
additively into a single parameter which we redefine again as a, for sake of simplicity; in effect we put "0 = 0. Then X 
is just the projected radius of the orbit which we denote by A. We therefore have just 3 search parameters for which 
a discrete mesh of filters is required: A = {k , A ,0, ,a) and A = (yl , i7 , a). The phase (4.2) becomes therefore: 



^ = KU + A cos{Q,u + a) 



(4.12) 



The metric gap, see Eq. (3.8) is now a 4 x 4 matrix, and we c ompu te it from Eq. (3.5) - with quantities scaled by the 
factor 2tt foT, and therefore replaced with 'tilde' - and Eq. ( 4.12 ). The first derivatives are simply: 



= cos{ilu + a) , 
$o = —Au sin(Slu + a) . 
I>a = — ylsin(riu + a) . 



(4.13) 



We need to compute the 10 time averages involving the above <tp {f3 = 0,1,2,3) and then compute the metric 
7jfc through the projection of gaf3 orthogonal to k, as discussed in Sec. [11. The exact analytical expression of 
Vp — det ||7jfe|| is very complex and not very illuminating. However, as discussed in the previous section, it is 
possible to compute det \ \jjk\\ in a closed form in the two relevant regimes: (i) fl ^ 1, the limit of several orbits 
during the observation time T, and (ii) SI ^ 1, the limit of monitoring a fraction of an orbit. 



1. The limit of large number of orbits 



We derive the e xpre ssion of the determinant of 7^ in the limit 3> 1. The elements of gajs and 7^^. are explicitly 
given in Appendix B 1. The proper volume element reads: 



\/det||7j-fc|| 



A^ 



96 



oiQ-\A^ 



(4.14) 



Fig. H shows the comparison of the analytical asymptotic expression ( 4.14 ) with the numerical evaluation of the full 
determinant. It is remarkable that for 10 - which corresponds to about 2 orbits completed during T - the two 
results are essentially identical. 



The proper volume Vp is given by integrating the volume element (4.14) within the appropriate limits. Here we 
need to integrate over A, n and a, within the ranges: 



4 ■ < A < A 

^miii — ^ — ^max ; 

< a < 27r . 



(4.15) 



Inserting Eq. (4.14) in Eq. (|3.3[), and taking into account the limits (4.15), we get the following proper scaled volume 
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Vv 



6%/6 



{A 



3 

max 



^min) (^11 



(4.16) 



Since the parameter space is 3-dimensional, = 3, to obtain the actual volume we need to multiply V-p by the factor 



6\/6 V 



(4.17) 



Notice that the factor 27rJ^ax/c is the maximum wavenumber of the gravitational wave that we want to detect. 
The main point in Eq. (4.17) is to observe how the volume scales. The number of filters increases linearly in the 
observation time T, so hierarchical searches based on T will not work effectively (compare this to the case of the all 
sky all frequency searches for isolated pulsars where the patches scale as T^). The volume is proportional to the cube 
of the size of the projected orbit along the line of sight Cp. So knowing, say, the inclination angle fairly well and the 
radius of the orbit will greatly reduce the computational load. Similar considerations apply to the frequency. 
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FIG. 2. Comparison of the asymptotic analytical expression and the numerical evaluation of the full determinant i n the large 
number of orbits approximation (f2 S> 1), for sources in circular orbits. The plot shows the proper volume element i^det ||7jfc|l, 
in arbitrary units, as a function of the dimensionless orbital frequency parameter fl = ujT. The solid line corresponds to the 
asymptotic expression in the regime f2 S> 1, Eq. (4.14), whereas the other curves correspond to the numerical evaluation of the 



full expression of det ||7jfc||, for several values of the initial phase of the orbit of the source: a = (dotted line), a = 7r/6 
(dashed line), a = vr/S (long-dashed line), and a — n/2 (dotted-dashed line). The size of the projected semi-major axis in 
dimensionless units is A — 10~^. Notice that for Q, > 20 the asymptotic expansion is indistinguishable from the full expression. 



In the case of blind searches, AA-' ~ ^innK^ ^O'' the parameters. Then the parameter volume is: 



V-p ~ 1.96 X 10^5 f J^) ( ^ 
V 1 kHz/ V 10" I 

and the total number of filters is given by: 



cm/ \ 10"*^ rad/sec 



J^ifi = 3%) ~ 1.25 X 10 
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Vv 

1015 



27r / 



T 



107 sec 



(4.18) 



(4.19) 
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This result indicates that the processing power to carry out a coherent bhnd search (even in the case of circular 
orbits) is outrageous; in fact, assume that the the costs are entirely dominated by the computation of the FFT and 
the power spectrum; then, each filter, over the entire search bandwidth, requires the following number of floating 
point operations: 



6 /in ax ^ 



l0g2(2/,„axT) + - 



(4.20) 



As a consequence the total processing power or computational speed S to keep up with the data is: 



S 



T 

2.1 X 10^*^ 



/max \ 

IkHzj 



Vio" 



MFlops . 



(4.21) 



We would like to stress that this estimation applies only to the orbital parameters and does not take into account the 
costs involved in searching for the source location in the sky and the spin down param eters. We would also like to 
point out the very steep dependence of the processing power, cfr. Eq. ( 4.18| ) and (4-21), on the maximum frequency 
/max upto which the search is conducted: S oc /^ax- This indicates that restricting the frequency band translates 
into a major saving in computational power. As one could have imagined, from previous results pertaining to isolated 
sources, the additional costs to search for the orbital parameters are absolutely prohibitive. This is also due to the 
fact, that the acceleration of a source with an orbital companion may be much higher than detector acceleration due to 
Earth's rotation and orbital motion, producing a much larger frequency drift which needs to be corrected for. Clearly, 
only prior information on the source parameters can make such a search feasible. In fact, when some parameters are 
known a priori, the number of templates to process the data reduces drastically. 
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FIG. 3. The errors on the orbital parameters required to carry out a search using only a few matched filters. The relative 
errors in the parameters Alj/lj (solid line), Aop/ap (dotted-dashed line ) and Aa/a (dashed line) are plotted for a constant 
thickness of = dl = 0.2 (A^ = 3 and fi = 0.03); the source corresponds to a NS/NS binary system in circular orbit; the 
parameters have been chosen to have the following values: lo — 1.027 x 10"* rad/sec, Op = 1.64 x 10^^ cm, a — 0.1 rad, and 

/max ~ 1 kHz. 
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It is interesting to estimate the thickness of the parameter space along the different directions as a function of the 
parameter values and errors, in order to gain insight into the accuracy that is required to carry out a search only with 
a few templates. For a 3 parameter search, and a mismatch jj, = 3% between the signal and the template, the filter 
separation, see Eq. (3.2), is dl — 0.2; it has to be compa red with the thickness T-', Eq. ( 3.18 ), for j = 1,2,3. Using 
our approximate expression of the metric, cfr. Appendix B 1, we obtain: 



T"^ ~ 1.49 X 10" 



Ac 



1km 



/n 



IkHz 



4.29 X 10" 



7.43 X 10" 



10-12 sec 

Aa 
10-5 rad 



10 



11 , 



( T 



( °P ) 
VlOiicm/ 



V 10^ sec 

/max \ 

IkHz J ■ 



/n 



IkHz 



(4.22) 



For typical binary neutron stars parameters, we also show in Fig. ^ the relative errors AA^/A-' that are required so 
that < dl. Our prior knowledge on the parameters must be very accurate in order to reduce the number of filters 
to only a few; typically an error AA-'/A^ < 10"^ (or smaller) is required. However, this is not at all uncommon for 
radio observations and this aspect will be discussed in Section ^ 



2. The fraction of one orbit case 

We consider now the opposite limit: f2 1. This corresponds to the physical situation where pulsar radio- 
astronomers usually apply the so-called " accelerated search" Q , as one monitors only a small fraction of the source 
orbital period, say less than one radian. The presence of a companion object for such a short time introduces an 
acceleration on the source motion that could be, as first approximation, treated as constant. 

For GW observations, the case ^ 1 is relevant for binary systems with very long orbital periods (of the order of 
one year or more) and T ~ 10^ sec, or when coherent integrations are performed on short data segments (say of the 
order of one hour). The latter situation is often encountered in hierarchical data analysis schemes, where one first 
uses matched filters on short time base-lines, and then concatenates together in a incoherent fashion the corrected 
data chunks. 

We derive an approximate expression for the determinant of the metric jjk by Taylor expanding it as a function of 
^ 1. The comput ation is trivial but cumb erso me, and we prese nt th e explicit expressions of the elements of the 
reduced metric, Eq. (|3.14 ), with phase model (4T) in the Appendix B2. We obtain the following volume element: 



det ||7,fc|| = (1 + ^"^^")^' ^2 ^8 ^ o(17^ A3) . (4.23) 



3628800V70 

It is interesting to note that for a = n/2 and 37r/2 - the NS is receding from, or approaching toward the detector 
essentially along the line of sight - the determinant tends to zero (at this order, but there are higher order corrections 
in f2 that depend on a, and do not vanish): in fact the signal appears Doppler-shifted by a constant indistinguishable 
off-set, depending on the velocity of the source, that one does not need to correct for. 



A comparison of the asymptotic limit (4.23) with the full expression of t he determinant is shown in Fig. M, where 



we present the results for two values of a — , 7r/3. Notice that Eq. ( 4.23 ) is very accurate upto 17 ~ 1 for a = 0, 
whereas it approximates correctly the full expression only for ;S 0.1 in the case a = 7r/3. This is due to the fact that 
when Of ^ 0, the determinant contains higher order correction terms in fi, that depend on a and vanish for a = 0. In 
the range 0.1 ;S ;S 1 they produce a sizable effect t hat o ne cannot ignore. 

By integrating Eq. (4.23) on the parameter range (4.15), and following steps similar to the previous case, we obtain 
the volume: 

= .....L^ (Ml^) ' [al^^^^ - a^^.„,J [(^„,..r)^ - (c..„,„r)^] . (4.24) 



24494400V35 

Following the same scheme as in the previous section, we can obtain the total volume of the parameter space, when 
no prior information about the source parameters is available: 

^~ ■ ^ 1 kHz y VlOiicm/ ^10-4rad/secy ^Tl^y 

^3.3x10- (J^^ (loT^) (lO-7rTd/sec) (l^^I^) ' ^"^'^^^ 
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The total number of filters, for a mismatch of 3% and the parameters values quoted in the previous expression is 
therefore N ~ 793 and 4.2 x 10^ respectively. The extra-computational burden to carry out searches of NS in binaries 
with long orbital periods as compared with the observation time, might be therefore affordable or even negligible in 
some cases. In fact, the processing power required to keep up with the data flow is respectively: 



5 ~ 1.4 X 10^ 



1.9 X 10" 



/n 



IkHz 



^) MFlops. 



(4.26) 




FIG. 4. Comparison of the asymptotic analytical expression and the numerical evaluation of the full determinant f or the case 
of observations lasting a fraction of the neutron-star orbit (fi <C 1). The plot shows the proper volume element det W'jjk \ \ , 
in arbitrary units, as a function of the angular orbital frequency, in u nits o f the observation time, = ljT. The solid lines 
correspond to the asymptotic expressions in the regime ^ 1, see Eq. (i.23), for selected values of the initial orbital phase of 
the source: a = and 7r/3 (see labels). For the same values of a we plot also -^/det ||7jfc| obtained by numerically evaluating 
the full expression (dotted line: a = 0; dashed-dotted line: a — vr/S). 



The costs are of course reduced if some orbital elements are known in advance, and the parameters need to be 
known only within a factor « 2 to allow a few templates search. In fact, the thickness of the parameter space in the 
three directions becomes now: 



l-96xlO- K§g|^ (ufe) [( 



7-6 X 10-3 (^) [(xep^) (li?)]' (^) 

Vl + cos(2«)] (,0^) (^) (lo^^)' (t^7)' (^) 

10-* soc-i ) ( IQii'^cm) ( 10-« soc-i ) (ihr) ( 1 kHz ) 



9.9 X 10-2 
2.5 X 10-2 



[a^(fc + i)7r] 
[a = (fc + i) tt] 
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- 3.8 X 10- (^) (^) [(tif^H^) (i^)] ' f 1^) [« = ^ -] (4.27) 

[6.3x10- (^) (^) [(to^^) (t^)]' (i^ml) 

where k is an integer. The difference with respect to the hmit 1 is striking; Fig. ^ clearly shows the steep 

dependence of the thickness of V on il, or, equivalently, the time of observation, where the behaviour of AX^ /X^ 
changes abruptly when T ^ P. Notice also the change of scaling of the T's depending on the location of the source 
on the orbit, i.e. the value of the parameter a. 



C. Elliptical orbits 



We have already described the orbit in Sections |IIB and IV A . As phase model we use Eq. ( 12 ) and take the 
expansion of the eccentric anomaly E' as a function of the mean anomaly M and the eccentricity e and truncate at 
the 7th order in e (see Eqs. ( |2.5| ) and Appendix This means we consider 7 harmonics in oj. The volume that 



we obtain is also correct upto this order in e. The derivatives of the phase (4.2) with respect to the six parameters 
K,,X,Y,e,uj,a, or, equivalently, k, A, e, cj, a, "0 are easily computed. We then proceed as in the circular orbit case, 
except now g^p is a 6 x 6 matrix and the 7^4, is a 5 x 5 matrix. The problem is therefore much more complex than 
in the circular case, and it is impossible to obtain a closed form expression of the determinant. However, we can still 
Taylor-expand the relevant expressions as done before, and keep the leading order terms. We present here a closed 
form expression of the volume element, and therefore the number of filters, in the asymptotic limit f2 3> 1. The limit 
r2 <C 1 can in principle be obtained by applying the same scheme adopted in the previous subsection. However the 
number of terms that we need to retain is so large that the computation becomes very cumbersome. We present only 
the numerical results in Fig. |5[ where the reader can derive the important scalings. 

In the limit f2 3> 1, using the parameters A, e, £7, a as coordinates, the volume element is given by. 



det I 



I7ife| 



^ r 33 

^ e- -e-^ 

32^/6 4 



41 
256^ 



cos 2-)/' 
32 



(4e3 - e^) 



cos 4-0 
256 ' 



(4.28) 



Comparing the former expression with the numerical evaluation of the full determinant (see Fig. it turns out - as 
it happened for e = 0, cfr. the previous sections - that Eq. ( 4.28D is accurate within a few percent even if the number 
of orbits completed during the time T is just about 2 or 3. 



In order to obtain the proper volume of the parameter space, we integrate Eq. (4.28) over the parameter range 



A ■ < A < A 

^min — ^ — ^max 7 
< P < P 

< a < 27r, 
< < 27r , 

and multiply by the scaling factor = {'^t^ fmaxT)^ . We have the following expression: 



Vp 



Fie, 



(Wn 



LOn 



where 



F(e) 



eHl 



41 
768* 



(4.29) 
(4.30) 



(4.31) 



(4.32) 



We can check how accurate this expression is by switching off the term in ( 4.32| ) and comparing with the fu ll one 



upto o(e°) . We find that the expressions agree very well to about e 0.7. From this we surmise that the Eqs. ( 4.31 ) 



and (4.32) will be correct upto e 
we found that up to e = 0.5 and e 



0.8. We also checked how well the leading order term approximates i^(e), and 
0.8 (as reference values), they agree within ~ 10% and ~ 35%, respectively. 
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FIG. 5. Comparison of the asymptotic analytical expression and the nume rical evalu ation of the proper volume element 
of the signal manifold for neutron stars in elliptical orbits. The plot shows det ||7jfc| |, in arbitrary units, as a function of 
the dimensionless orbital frequency parameter Q = uiT. The solid line corresponds to the exact expression of t he determinant, 
computed numerically. The dotted line refers to the asymptotic expression of the determinant, cfr. Eq. ( |4.2^ ), in the regime 

S> 1; the dotted-dashed line is the asymptotic expression of the determinant for f2 ^ 1, in which case det | \yjk\\ oc Q^^ . 
The size of the projected semi-major axis in dimensionless units is j4 = lO"'^; the other orbital elements are chosen as e = 0.5, 
and a = if) — 0. 



For blind searches, the size of the parameter space is: 



where, we have evaluated Eq. (4.31) for e^ax = 0.5. For Cmax = 0.1 and emax — 0.8 the volume changes by the factors 



0.04 and ~ 2, respectively. The total number of filters now is: 



AA(/i = 3%) ~ 1.1 X 10^6 ( ) . (4.34) 



,10' 

Indeed, for e 7^ 0, prior information about the source parameters are essential in order to reduce the computational 



costs, see Eq. (4.21). In particular the steep dependence on Op and /, both to the fifth power, indicates that the 
computational load could be cut down substantially if we knew them fairly accurately. These are the crucial parameters 
which govern the computational costs. Af still scales linearly with ujT. For sufficiently small values of the eccentricity, 
it scales quadratical ly wi th e. 



As we did in Sec. tVB , we can estimate the thickness of the parameter sp ace along a given direction and compare 
it to the filter separation dl, which for A'' = 5 and fj, = 3% is ~ 0.15 (see Eq. (^)). This is useful in obtaining insights 
on the accuracy needed on the parameters for reducing the computational costs by a large factor, so that they become 
affordable. In Fig. ^ we show the relative errors AA-'/A-' , for typical orbital elements of a NS/NS binary, as a function 
of the observation time, where we set = dl. For observation times involving several orbits, 0^1, one needs to 
know a parameter with an accuracy of the order of one part in a million, in order that the search is carried out with 
the help of only a few filters. However, as in the circular orbit case, AA-'/A-' increases very steeply when T becomes of 
the order of P, or shorter. In this case a 1% accuracy (compare it to the case e = 0, where AA-'/A-' ~ 2) is sufficient 
to carry out the search with few filters. 
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FIG. 6. The errors on the orbital parameters required to carry out a search using only a few matched filters for a binary in 
eccentric orbit. The relative errors in the parameters Aoj/oj (solid line), /S.X/X (dotted-dashed line), Ae/e (dotted line) and 
l\a/a (dashed line) are plotted for constant thickness equal to the filter separation dl ~ 0.15, corresponding to = 0.03, for 
an edge-on eccentric orbit with parameters: u — 1.027 x 10~^ rad/sec, Op = 1.64 x 10^^ cm, e = 0.5. 



V. COMPUTATIONAL COSTS FOR KNOWN NEUTRON STARS 



In the previous sections we have analyzed the extra computational costs required in correcting for the Doppler effect 
produced by the motion of a NS around a companion object, on a monochromatic GW signal. We are therefore in a 
position to apply the former results to known neutron stars which are amenable to searches by Earth-based detectors. 
Here we deal with two main classes of sources: the NS's known as radio pulsars and the NS's in LMXB's, in particular 
Sco X-1. 



A. Radio pulsars 

The last published version of the pulsar catalogue contains 706 NS's, of which 44 are in binary systems. 
Actually, the total number of pulsars known today is ~ 1200, thanks mainly to the recent multi-beam survey of the 
southern sky carried out at the Parkes radio-telescope in New South Wales, Australia p^-|30[|. A further considerable 
increase in the number of pulsars is expected over the next few years when the upgrade of other radio telescopes will 
be completed, allowing a survey of the northern sky similar to the one carried out at Parkes. However, the newly 
discovered sources have not yet been included into a publicly available data base, and we therefore restrict ourselves to 
the catalogue of 1995. Clearly, when the newly discovered NS's, and their fundamental parameters will be available, 
our present analysis can be applied to the new objects in a straightforward way 

In order to compute the number of filters required to target the 44 NS's, we apply the scheme presented in the 
previous sections. More specifically, we use the full and exact expression of the parameter space volume computed 
numerically, with the values of the parameters and errors quoted in [TtI ; the fundamental ones are also shown in 
Table |. For several sources the orbital parameters are known very accurately, so that we find that the search is 
limited either over a small region of the parameter space, or to a smaller number of dimensions (< 5), or both. We 



first estimate the actual number of dimensions of the signal manifold, following the discussion in Section IIIB, and 
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then calculate the total number of filters from Eg. (3.27). We set the value oi 5 = 0.1 - the thickness of the parameter 
space in units of the filter distance, see Eq. ( 3. 22 ) - which distinguishes essentially the known parameters from the 
questionable ones. 



TABLE III. The number of filters required to search for NS's known as radio pulsars in binaries. We list the total number 
of filters Af required to search for the 44 NS's in binary system contained in the Taylor et al. catalogue for an observation 
time T = 1 yr and a maximum mismatch of 3%. 



PSR J 


PSR B 


Af 


0024-7204E 


0021-72E 


1.8 X 10^ 


0024-72041 


0021-721 


4.1 X lO'^ 


0024-7204J 


0021-72J 


9.2 X 10^ 


0034-0534 




1. 


0045-7319 




signal outside the sensitivity window 


0218-^4232 




1.5 X 10^ 


0437-4715 




1 


0613-0200 




1 


0700+6418 


0655+64 


1 


0751-^1807 




1.1 X lO'' 


0823-^0159 


0820+02 


signal outside the sensitivity window 


1012-1-5307 




4 


1022+10 




2 


1045-4509 




16 


1300+1240 


1257+12 


1 


1302-6350 


1259-63 


2 


1312+1810 


1310+18 


3.9 X 10^ 


1455-3330 




2.6 


1518+4904 




2 


1537+1155 


1534+12 


1 


1623-2631 


1620-26 


8 


1640+2224 




4 


1641+3627B 


1639+36B 


6.4 X 10^ 


1643-1224 




7.3 X 10* 


1713+0747 




1 


1721-1936 


1718-19 


signal outside the sensitivity window 


1748-2446A 


1744-24A 


4 


1803-2712 


1800-27 


signal outside the sensitivity window 


1804-0735 


1802-07 


4 


1804-2718 




1.6 X 10'' 


1823-1115 


1820-11 


signal outside the sensitivity window 


1834-0010 


1831-00 


signal outside the sensitivity window 


ioO ( +uy4o 


1 Q 1 An 

iooo+uy 


i 


1910+0004 




1.2 X 10^ 


1915+1606 


1913+16 


1 


1955+2908 


1953+29 


4 


1959+2048 


1957+20 


2 


2019+2425 




1 


2033+17 




8.9 X 10" 


2130+1210C 


2127+llC 


33 


2145-0750 




4 


2229+2643 




16 


2305+4707 


2303+46 


signal outside the sensitivity window 


2317+1439 




1 



The results of our analysis are presented in Table III , which lists the number of templates required for an integration 
time of T = 1 yr. The key point to notice is that the most interesting NS's require at most a few filters in correcting for 
the orbital motion: J0437-4715: NS/white dwarf system; J0437-4715: NS/white dwarf system; B1257+12: NS/planets 
system; J1537+1155: NS/NS binary; J1623-2631: NS/white dwarf system; and J1300+1240: NS/main sequence star 
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binary. The upper-limits on he are, in fact, within the sensitivity of the first or second generation of detectors, cfr. 
Figure |l]. The data analysis costs for such sources is therefore minimal. Also considering that the frequency is very 
accurately known within a small bandwidth, one can substantially down-sample the data stream around the expected 
source frequency. Moreover, the orbital companions of the above mentioned five neutron stars are also very different, 
which boosts the astrophysical interest in monitoring them. 

On the other hand, for those NS's about which no information about the spin-down is available, we cannot set 
upper-limits on the signal strain and hence the computational costs become considerable. The number of filters is 
in the range ~ 10^ — 10^", for a year's worth of integration time. Therefore, even assuming a search in a narrow 
frequency range, around twice the radio frequency, the computational burden is overwhelming. The situation is made 
worse by the lack of information about / and higher time derivatives of the frequency, so that one would also have to 
scan the space of spin-down parameters, increasing the costs further. It is clear that more accurate data from radio 
observations are necessary in order to narrow down the uncertainty in the source parameters for making such searches 
feasible. 

To summarize, the interesting radio pulsars in binary systems do not lead to large extra computational load when 
correcting for the orbital parameters of the source. On the other hand, for the NS's for which we have no idea about 
the GW strain amplitude, the computational costs are not affordable. Such sources call for a follow-up of radio 
observations for obtaining more accurate measurements of the orbital parameters. 



B. Sco X-1 



Low Mass X-ray Binaries (LMXB's) could be primary GW sources for Earth-based detectors (see discussion in 
the Introduction). Our current astrophysical understanding and signal modelling suggests that at least one system, 
namely, Sco X-1 - the brightest X-ray source in the sky, located at a distance 2.8 ±0.3 kpc - is just within the reach 
of the first generation of detectors with very long coherent observation times (of the order of 2 years), and clearly 
detectable by the second generation of instruments (in particular LIGO II, in narrow-band configuration), for an 
integration time as little as « 15 days. 

Sco X-1 orbits a low-mass companion in an essentially circular orbit with a period P = 0.787313 (1) days; the 
orbital angular velocity is a; ~ 9.24 x 10~^ rad/sec and is known with an error Auj ~ 2.34 x 10"^*^ rad/sec. The 
radial component of the orbital velocity is = 58.2(3.0) km/sec therefore the projection of the semi-major 
axis along the line of sight is Op = Wr/w ~ 6.3 x 10^" cm. The relative error on the semi-major axis is Aop/cp = 
AP/P + Av-c/v-c ~ Av-c/v-c, since the period is very accurately measured. This yields an error Aa^ ~ 6.5 x 10^ cm. 
The position of the NS in the orbit is known with an error Aa = 0.2 rad p3[ . The GW emission frequency is not 
very well determined; conservatively we assume it to be within the band 500 Hz - 600 Hz, and therefore we take the 
bandwidth A/ « 100 Hz In the expressions for thicknesses of the parameter space below, it is appropriate to 
replace /max by the bandwidth A/ over which the search would be launched. 

The main problems that one faces in detecting (or setting upper limits on) GW's emitted by Sco X-1 are: (i) the 
drift in the GW frequency, which is not likely to be monotonic, and is related to the time-varying accretion rate 
which is ultimately responsible for the net quadrupole moment, and (ii) the need for taking into account the Doppler 
phase modulation caused by the orbital motion. Present order-of-magnitude estimates of the frequency evolution 
suggest that for a time duration of w 10 days the signal is confined within a single frequency bin; it determines 
the longest coherent integration time (or base-line duration of observations) that one can employ in searching for a 
monochromatic signal. As a conservative limit, we take a week: T — 6.05 x 10^ sec. Notice that during this time the 
radiation is monochromatic, but emitted in an unknown frequency bin inside the band 500 Hz ~ 600 Hz; the position 
of the source is also known very accurately, so that one is essentially left with the orbital parameters that one must 
search over. The estimate of the computational costs that we present therefore reflects the total costs involved in 
searching for Sco X-1. 

It is interesting to consider first the thickness of the parameter space and the conditions such that Sco X-1 is a 
one-filter target. Notice first that = 1 for T — 1.08 x lO^sec — 3.01 hrs. In the limit ^ 1, therefore up to ~ 3 
hours, we have: 



~ 9.3 X 10" 



T'^ ~ 1.6 X 10" 



yjl + cos(2a) 
5 — cos(2q;) 



Vl + cos(2a) 



Aa„ 



, 6.5 X 10^ cm 
Aw 
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9.24 X 10-5sec"i / V Ih^ 



A/ 
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Notice that the expressions for T^p and are given for a ^ fc + 7r/2, where k is an integer; the expression for T" 
is valid if a ^ fc7r/2. In t he ot her cases one can easily derive the value of the thickness by applying the appropriate 
formula, as given in Eq. ( 4.27 ). 

If one takes an hour long data stretch, Sco X-1 in a one-filter target. For longer observational times up to T « 3 
hours, the search requires at most a few templates to take care of "border effects" . A coherent search over a bandwidth 
~ 100 Hz can be performed very inexp ensiv ely, and one needs a fraction of an MFlop to keep up with the data flow, 
which can be easily derived from Eq. (4.21) with /max replaced by A/ « 100 Hz. However, the computational cost 
increases steeply (as shown in the previous section) when T becomes of the order of the orbital period or longer. In 
fact, for ^ ^ 1, the thickness of the parameter space becomes: 



~ 96. 



~ 6.3 X 10~3 



6.5 X 109 cm/ V 100 Hz 



T" ~ 93.6 



Ac 



2.3 X 10-10 sec-i / V 6.3x101° cm 



A/ 



T 



10^ sec 



A/ 
100 Hz 



0.2 rad/ V 6.3x101° cm / V 100 Hz 



(5.2) 



One needs a mesh of filters for ap and a, but up to T w 1 month the orbital period is not a search parameter (one 
might actually need up to 2 filters to take care of the error s in P ). The reduced metric 7jfc which we need to consider 
is given by a 2 x 2 matrix on Op and a. By applying Eq. ( ^.25| ), with /i ~ 0.03 and T ^ 1 month, we obtain. 



J^(N= 



iN=2) 



5.2 X 10*= 



(5.3) 



A one week coherent search involving 5 x 10^ filters requires « 84 GFlops of processing power to keep up with the 
data flow. For longer integration times T ^ 1 month, uj becomes a search parameter; notice that in this case the 
signal is not anymore monochromatic, and one needs to search also over the spin-downs parameter space (wh ich is 
not trivial, due to our poor understanding of the accretion rate and angular momentum transfer, see Sec. IB); the 
number of filters to correct only for the Doppler shift induced by Sco X-1 orbital motion is given now by: 



J^(N= 



(N=3) 



3.1 X 10" 



T 



lO^sec 



(5.4) 



By setting a coarser grid, increasing the mismatch from, say, 3% to 30%, will reduce the number of filters in the two 
cases (5.3) and (5.4) by a factor ~ 10 and ~ 30, respectively. The values still indicate that the computational costs are 
too high which is also due to the large frequency band that one needs to search through. These results clearly show 
that it is mandatory to determine the orbital parameters more accurately in order to reduce the parameter volume 
so that Sco X-1 can be monitored continuously. Eq. (5.2) implies that the errors on Op and a must be reduced by 
3 orders-of-magnitude, so that only a few templates are required for the search. Although this calls for a very large 
observational step, it is still possible. 



VI. CONCLUSIONS 



We have estimated the extra-computational costs required to search for CW sources in binary systems. This work 
was motivated by the fact that LMXB's, a key class of sources, possibly, for the initial generation of detectors and 
certainly for the enhanced detectors, are NS's orbiting a binary companion. Also several NS's detected in the radio 
band are in binaries and will be continously monitored by the Earth-based network of GW detectors. In order to 
disentangle the computational costs arising from correcting for the Doppler phase modulation produced by the source 
motion around the orbital companion, we have assumed that the signal is exactly monochromatic and the location of 
the source in the sky is perfectly known. Our results provide an estimate of the computational burden associated with 
these extra dimensions of the parameter space which has been neglected so far (for rather obvious reasons). We have 
deduced closed form analytical expressions for the number of templates required to carry out the analysis. We have 
obtained the scaling of the computational costs as a function of the orbital elements and the observational parameters 
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(the source emission frequency and the observation time); they are ready-to-use tools to evaluate the computational 
resources needed and the trade off for search strategies. We have applied our analysis to known radio pulsars and Sco 
X-1, estimating the computational costs involved in targeting these systems. The analysis also addresses the need of 
further astronomical observations of known NS's, with the aim of reducing the errors in the source parameters, so that 
the search is made inexpensive. In fact, a byproduct of our analysis is the rigorous formulation, in the geometrical 
framework of data analysis, of the conditions, i.e. the size of the parameter space, such that a NS can be considered 
a one-filter target; our result is completely general, and can be applied to any class of sources (not necessarily CW 
sources) whose waveforms can be modelled reliably. 

Our analysis is clearly limited in several respects. Possibly, the main accomplishment is a quantitative understanding 
of the key issues that require intensive attack, and identifying the bottlenecks that prevent us from carrying out blind 
searches of NS in binary systems. There are three clear directions of future work that we believe should be pursued 
vigorously, particularly in the light of the current data analysis effort in setting up GW search codes: 

1. The software implementation of matched filter based searches of a NS in a binary system, for given values of the 
source parameters such as location in the sky, spin-down(s), orbital elements. This task is simple in principle, 
as the fundamental building blocks are already available: search codes for GW's from isolated NS's and publicly 
available software (such as TEMPO) in the radio pulsar community to correct for the Doppler effect of the signal 
phase for a source in a binary system. However, the implementation of this might not prove to be so trivial. 

2. The estimation of the computational costs keeping into account both the orbital elements and the the source 
position and the spin-down parameters. This is the composite problem where the full expression for the total 
costs needs to be obtained. It will take into account possible correlations between the search parameters 
investigated here, and those that describe isolated neutron stars. 

3. The investigation of suitable hierarchical strategies, that can reduce the computational costs as compared to 
fully coherent searches. The current vigorous effort in the study, and software implementation, of hierarchical 
algorithms for isolated sources, the so-called stack-and-slide and Hough-pattern-tracking searches provide already 
a template for investigations for a larger parameter space. Moreover, radio and X-ray astronomers are currently 
investigating algorithms suitable to carry out analogous searches in the electro-magnetic spectrum [^,^ . 

One limitation of our approach is that we have completely neglected relativistic effects in the model of a binary 
orbit. In particular, the advance of the periapse is very significant in relativistic systems (the largest periapse advance 
reported in p7[ | is ~ 4 deg/yr). This, and other complications, with respect to the simple Keplerian model that we 
have adopted here, need to be taken into account in future studies, and addressed carefully. 

We would like to conclude by stressing that the data analysis problem that we have addressed here is common to 
electro-magnetic observations of pulsating stars in binary orbits, where typical examples are radio and X-rays binaries. 
We can therefore hope that the current efforts in other related research fields can help in providing a breakthrough 
in this area. 
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APPENDIX A: COEFFICIENTS OF THE TAYLOR EXPANSION OF THE ECCENTRIC ANOMALY 



Here we give the expressions of the coefficients Cfc(e) and Sk{e) which appear in the expressions for cosE and 
Vl - e2 sinE in Eq. (pj]), upto the 7-th order in e [|||: 
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(A2) 



we refer the reader to p3|, or any other standard book of celestial mechanics for the general expressions of the 
coefficients Cfe(e) and Sk(e). 

APPENDIX B: THE APPROXIMATE METRIC COMPONENTS FOR CIRCULAR ORBITS 



We pr ovid e t he t he exp ressio ns, to the relevant order in A and fl, of the components of th e metrics gap and "fjk, 
see Eqs. (3.8), ( |3.9D , (3.10) and (3.12), for a NS in circular orbit, using the phase model (4.12) in the two asymptotic 
limits SI > 1 and SI < 1. 



1. The asymptotic limit of several orbits 



We derive the expressions of the components of ga/s {a, P — 0,1, 2, 3), and 'jjk («, j — 1, 2, 3) given by Eq. ( 3.12| ) for 
S7 ^ 1. We recall that, in general, A <^ The time derivatives are given in Eq. ( 4.13| ), and the time averages are 
computed through Eq. (p.9|). 

To the leading order in A and the elements of gap read: 
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Substituting these components into Eq. (3.10), which is related to jjk through Eq. (3.12), one derives the foUowing 
expressions, to the leading order in A and 17: 
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(B12) 
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(B14) 
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(B16) 



The determinant (4.14) is obtained from the latter components by retaining the leading order terms in A and 1/0. 



2. The asymptotic limit of a fraction of one orbit 

The expression of the determinant in the case 0^1 can be simply derived by Taylor expanding the components 
of ga0 and jjk in powers of O. Here we present the expressions of gap upto 0^°: 
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Substituting the former components into Eq. (3.1C), which is related to ^jk through Eq. ( ^3.12 ), one derives the 
foUowing expressions for the components of the metric jjk upto order fl^'^: 
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The determinant (4.22) is then derived from the latter components and Taylor expanding the resulting expression to 
the leading order in Q. 



APPENDIX C: THE APPROXIMATE METRIC COMPONENTS FOR ELLIPTIC ORBITS 

Below we list the metric components cjajs for a > /3, in the coordinates k, X,Y,e,ip,cx upto order and in the 
leading order X, Y for the case 1: 
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In the previous expressions we have also dropped all the terms of order and higher, and set a — Q. Substituting 
these components gjk into Eq. ( 3.10| ), and consistently working withi n the appropriate approximations, one can derive 
(to the leading order in e^), the expression for the determinant in (4.28). 
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